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We theoretically investigate a quantum dot coupled to fermionic (electronic) leads and show how 
zero-point quantum fluctuations stemming from bosonic environments permit the rectification of 
the current. The bosonic baths are either external impedances modeled as tunable transmission 
lines or LC resonators (single-mode cavities). Voltage fluctuations stemming from the external 
impedances at zero temperature are described through harmonic oscillators (photon-like excitations) 
then producing the quantum vacuum fluctuations. The differing sizes of the zero-point fluctuations of 
the quantum vacuum break the spatial symmetry of the system if the quantum dot is coupled to two 
reservoirs or two junctions with different bosonic environments. We consider current rectification 
and power production when the system is operated as a heat engine in both non-resonant and 
resonant sequential tunneling cases. 


I. INTRODUCTION 

The dream of using fluctuations from the quantum vac¬ 
uum as a limitless power source has inspired many, from 
science fiction writers, to scientists. While it is certainly 
possible to extract energy from the vacuum (a pair of 
Casimir plates moving toward each other can do exter¬ 
nal worlP), the difficult part is to find a way to do it 
cyclically or continuously without putting in more energy 
than you get out. Indeed, if we consider the quantum 
vacuum as the many body ground state, then it is clear 
that if the Hamiltonian itself is not changed, no energy 
can come from it. Quantum vacuum fluctuations arise 
simply because the operators in question (position, mo¬ 
mentum, particle number, etc.) do not commute with the 
many body Hamiltonian, and are no longer good quan¬ 
tum numbers. 

We consider the possibility of rectification of electri¬ 
cal current from quantum vacuum fluctuations for the 
specific case of electrical transport through a quantum 
dot, and its use as a nano-heat engine. The subject 
of quantum thermodynamics has been rapidly growing 
in interest! 2 . This has been driven both by fundamen¬ 
tal questions as well as possible applications for energy 
harvesting and cooling. There has been much recent ac¬ 
tivity on thermoelectrics and quantum dots 3 9 and also 
pioneering experimental works 10 13 . 

In this work, we consider the systems displayed in 
Fig. l(a,b) and show how the zero-point quantum fluc¬ 
tuations can break the symmetry between left and right 
then producing a direct current without the application 
of a bias voltage across the system. The left and right 
leads contain electrons, tunnel coupled to one another 
via a quantum dot. The dot is described with coupling 
capacitances Cp, Cr. Additionally, this system is part of 
an electrical circuit with finite impedances Zp and Zr 
connected to the system via the capacitances C£ and C£ 
respectively (and an additional impedance Zy coupled to 
the dot through the capcitance Cy in the case of the setup 
of Fig. la). The impedances can be modeled by a chain 
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FIG. 1. (color online). Quantum vacuum nano-engines with a 
quantum dot in the Coulomb blockade regime. The panels (a) 
and (b) represent the circuits of the two distinct setups. The 
panels (c) and (d) are pictorial representations of the tun¬ 
neling processes respectively for the setup (a) and (b): the 
lateral parts correspond to the energy diagram of the elec¬ 
trons in the corresponding lead (at zero temperature for (c) 
and at non-zero temperature for (d)) and the central part 
represents the energy levels of the quantum dot. The straight 
dotted arrows represent a tunneling event from left to right, 
associated to curvy arrows which correspond to photon emis¬ 
sion/absorption. (a) and (c) Photons from a hot impedance 
Zd excite zero-temperature electrons tunneling onto the quan¬ 
tum dot. These electrons emit photons in the cold impedances 
to tunnel off of the quantum dot. (b) and (d) Electrons at 
finite temperature tunnel through the quantum dot and emit 
photons in the cold impedances to tunnel on and off of the 
quantum dot. In both cases, the asymmetry in the quan¬ 
tum vacuum fluctuations from the zero temperature latteral 
impedances leads to rectified current. 

of LC oscillators, which when quantized, behave as a sys¬ 
tem of quantum harmonic oscillators, whose frequency c Oi 
depends on the impedance and inductance of the chain. 
The vacuum energy of the chain is Eo = (1/2) JA hooi, 
and leads to fluctuations of the charge (or voltage) on the 
plate of the coupling capacitor to the quantum dot 14 . 
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We will address resonant and non-resonant sequential 
tunneling processes, associated to the emission of pho¬ 
tons in the zero-temperature transmission lines, and we 
assume spin-polarized electrons for simplicity. For the 
regimes we consider, the spin just multiplies by two the 
number of conducting channels. The distinct regime of 
large impedances (dynamical Coulomb blockade) was ad¬ 
dressed in Ref. m We will also introduce a different 
setup, where the leads are now coupled to frequency 
selective cavities (or tunable LC reonators which can 
be built with current technolog^ 16 -!). These two classes 
of bosonic environments coupling to mesoscopic circuits 
have attracted some attention both experimentally 17 27 ^ 
and theoretical!}®®] Exotic many-body physics has 
also been addressed^®! 

We now briefly discuss the physics of our setup. 
Naturally, no current is rectified when the system is put 
in state of overall zero temperature: it is necessary to 
give some energy to the system for a current to flow. We 
consider two ways to do this: first, we consider the case 
when the conduction electrons are at zero temperature 
so that the charging energy of the quantum dot forbids 
electrons to tunnel onto it. However, if the quantum dot 
is capacitively coupled to an electrical environment via 
capacitance Cd with impedance Z^, and temperature T, 
the electrons have a new way of gaining energy: photons 
from the hot impedance can be absorbed, allowing the 
electrons to tunnel onto the dot. This corresponds to 
the setup displayed in the Fig. la. Second, we consider a 
model where the electrons in the two leads are in thermal 
equilibrium, both at the same small temperature T. 
This corresponds to the setup displayed in the Fig. lb. 
In both cases, the finite charging energy of the quantum 
dot requires an energy E from the incoming electron 
to be occupied. In the first case, this is supplied by 
absorbing a photon from the hot circuit; in the second it 
is from the fact that some of the electrons are thermally 
populated into that energy. Giving some energy to the 
system is however not sufficient to induce a current. We 
also need to break the spatial symmetry of the system 
in order to favor one direction of tunneling events with 
respect to the other. This can be done by coupling 
the leads to zero temperature electrical circuits with 
different impedances Zl and Zr (or the same impedance 
but with different coupling capacitances Cl and Cr). 
To leave the dot, the electron must reenter the Fermi 
sea of the quasi-particles surrounding it. This process 
is sensitive to the vacuum fluctuations of the electrical 
circuit. If the impedances coupled to the two leads 
are different, the fluctuations of the vacuum will break 
Left/Right symmetry, giving a preferred direction for 
the particle to exit the system into - thus giving current 
rectification. 

The paper is organized as follows. In Sec. II, we intro¬ 
duce the quantum many-body physics corresponding to 
Fig. 1, and discuss the physics in terms of P(E) theory®. 
In Sec. Ill, we derive the quantum vacuum-induced rec¬ 


tified current for the different geometries considered, and 
show its dependence on the physical parameters of the 
system. In Sec. IV, we add a load voltage to the sys¬ 
tem, and consider its operation as a heat engine, finding 
the power produced and thermodynamic efficiency versus 
temperature and asymmetry. We give our conclusions in 
Sec. V. In Appendix A, we quantize the voltage fluctua¬ 
tions through a transmission line description. 

II. ANALYSIS 

We let Ck be a fermionic annihilation operator on the 
left lead of momentum fc, similarly, d and c q describe 
electrons on the dot, and the right lead. 

The system Hamiltonian is, H = Hr + Hr + Hd + 
H ph + H c , where H L = Ld e k + e5V L )c\c k , Hr = 
+ eSVR)c^ q c q are the Hamiltonians for the left and 
right leads with electron energies e/~ and e q . H p l is the 
Hamiltonian of the environment in isolation and SVl,r 
are the voltage fluctuations on the capacitors induced by 
the presence of the Left/Right environments. A detailled 
description of these terms is provided in the appendix 
A. We define the zero of energy as the Fermi level of 
the leads. The quantum dot Hamiltonian, describing a 
single quantum level, is Hd = Ed + E c (d) d — 1). This 
specific form of charging energy assumes spin-polarized 
electrons. Ed corresponds to the energy level of the dot, 
as defined in Fig. 1 and E c > 0 is the charging en¬ 
ergy of this level, which corresponds to the differe nce of 
energy when the dot is charged, versus uncharge d 45 1 46 1 . 
The external impedances are coupled to the electrical 
circuit through coupling capacitances, allowing then an 
efficient thermal isolation. We find by a circuit analysis 
that E c = e 2 /2Cs — eJ^jiCj/C^SVj, where the sum is 
over L,R, d. Here C e = Cl + Cr + Cd and SVd cor¬ 
responds to the voltage fluctuations linked to the cen¬ 
tral impedance d, which is only present in the case of 
Fig. 1(a). The only effect of these impedances is then to 
induce environment-assisted tunneling events. 

The tunneling terms are He = ^2 k tkdd^Ck + h.c. + 
E qtdqCqd + h.c.. Here, tkd is the tunneling matrix el¬ 
ement onto the dot from the left lead with momentum 
&, td q is the tunneling matrix element from the dot to 
the right lead with momentum q. The effect of the fluc¬ 
tuating voltages on the tunnel junctions and coupling 
capacitors can be taken into account by a unitary trans- 
formation U = exp (i<j> L ]T fc c\c k +i(f> R J2 q c\c q +i<j) d d)d). 
This transformation has the effect of eliminating the fluc¬ 
tuating potentials on the bare system Hamiltonians, and 
putting a fluctuating phase on the tunneling elements, 

tkd -t t k d exp (i$ L ),td q -t t dq exp(-i$i{). (1) 

Here, we have introduced the total phase = <fii — <fid, 
where pi = (e/ti) dt'6Vi(t') for l G {L,R}, and 
Pd = (e/d) /Cz) f* dt'SVj(t') where the sum is 

over j G {L,R, d}. The voltage fluctuations SVj are 
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linked to the environmental degrees of freedom whose dy¬ 
namics is governed by H env , as shown in the appendix. 
Below we show how the impedances in the different cir¬ 
cuits of Fig. la and Fig. lb contribute to the noise. 

We treat the problem within P(E) theory®, assuming 
that tunneling on and off of the dot can be described 
with tunneling rates described within a golden rule pic¬ 
ture in perturbation theory. This also assumes that envi¬ 
ronments relax faster than the tunneling time scale rates. 
In this picture, we do not take into account the backac- 
tion of the electrons on the excitations of the impedances. 
The tunneling rates are controlled by the probability of 
the electrical circuit to either absorb or emit a photon. 
This probability for lead l E {L, R} is given by 

/ oo 

dte iEt / h e Kl (*), ( 2 ) 

-OO 

where Ki(t ) = — <£j(0))<^(0)) (where the 

variables are defined above). 


In the case of cold electrons and hot photons (Fig. 1(a)) 
we have K t = (r$ + rf m )ki + vjkd for (l,m) = 

(L, R), (R, L) where we have introduced the coupling 
constants to the environmental baths, rju = ( C m + 
Cd)/Cx,rjim = Cm/ Ce, r] d = Cd/C E (m = R when 
l = L and vice versa). ki(t) and kd(t) are the lead 
and dot correlation functions taken in isolation kj(t) = 
((<; f)j(t ) — 0 j(O)) 0 j(O)) and they read 
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where Zi = yjLi/Cf is the impedance coupled to the 
lead l. Zi is expressed as a function of Li and Cf, which 
are the values of the inductances and capacitances which 
compose the transmission line (see appendix A). Physi¬ 
cally, the correspondence between impedance and trans¬ 
mission line just fixes the ratio between inductances and 
capacitances, therefore we can choose the capacitances 
to be equal to Cf and adjust the inductances such that 
the ratio y/Li/Cf gives the physical resistance at uo = 0 . 
R q = h/e 2 is the resistance quantum. 

For the sake of simplicity in this case of the setup 
1 (a), we have specified the case for which both left 
and right leads have the same correlation function 
ki,(t) = &r(£) = k 0 (t) in isolation at zero temperature, 
and the asymmetry of the system comes from different 
values of Cl and Cr. This is one particular choice 
to break the Left/Right symmetry, which leads to the 
simple expression above for Ki. There are other ways to 
break this symmetry such as having different impedances 
Zl and Zr. The important point is that the capacitive 
coupling to the hot impedance now enters the Pi(E) 


functions for the tunneling electrons. 

In the case of hot electrons and cold photons (Fig. 1(b)) 
we have rjd = 0 and for Cl = Cr we find that Ki = 
ki given directly by the first line of Eq. <[3]). In this 
case we will take distinct impedances leading to different 
functions Kr^) and ATr(£). 


III. RESULTS FOR THE DIFFERENT 
GEOMETRIES 

A. Cold electrons and hot impedance Zd (Fig. 1(a)) 

The effects of ko(t) and kd(t) are quite different and 
we will first consider their effects independently. As 
stated previsouly we consider for ko(t) the case where 
the zero frequency external impedance of the T = 0 
transmission lines Z = R = \JL/C c is small compared 
to the resistance quantum, R R q = h/e 2 , and de¬ 
fine the large a = R q /R 1 (where L = Lr = Lr 
and C c = C£ = Cr). In this case, we have for 
long time, ko(t) = —(2/a)ln(aE c t/7rh) + iir/2 + q e , 
where q e is the Euler constant 44 ! The isolated ef¬ 
fect of this cold impedance would lead to a Pq(E) = 
(1/ft-) f^° dt exp (iEt/ti) exp [ko(t)] function that van¬ 
ishes for negative energies, and has a power-law diver¬ 
gence for small positive energies, signaling the onset of 
the orthogonality catastrophe. In contrast for the hot 
impedance, still in the small resistance limit, the cor¬ 
relation function becomes kd(t) = — ( ReZd/R q )kRTt/h 
in the long time limit. The isolated effect of this cold 
impedance would lead to a normalized Lorentzian for 
Pd(E) = (1/ft) dt exp (iEt/h) exp [kd(t)] of width 
A = (ReZ d /R q )k B TW 

The combined Pi(E) functions from both effects is 
found by taking a convolution of the two Po(E) and 
Pd{E) functions (convolution theorem), provided the 
coupling constants rji to the various circuits are prop¬ 
erly accounted for. This turns out to be straightforward 
because they can be absorbed into effective impedances 
Oil = ot/(r)fi + Vim), Charging energies E C>1 = (Vu+V? m ) E c 
and Lorentzian widths, = A 77 ^, giving 

yp2/o;; g 2 / y e /ai ~ j2 . 

PAE) = -~--—Im(— iAi— E) & i 

U ' Y(2/a l )^(2^/a l )(a l E c ^/^ V 7 

(4) 

The Left/Right asymmetry is now solely contained in 
the effective parameters fti/ cxr. These Pi(E) functions 
permit us to calculate the tunneling rates T±j between 
the leads and the dot, from left to right (+) and right to 
left (-), 


T±, L m = t 0 ,l j 

! def(e±nT^/2)P L (e), 

(5) 

= Rut 

f def(eT^T ld/2)P R (e). 

( 6 ) 
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Here T 0? / = 2ii\ti\ 2 viVi is the bare tunneling rate given 
in terms of the tunnel coupling £/, the density of states vi 
and the volume Vi of the lead, and we have included the 
possibility of a bias V (/x = eV) on the system. Q de¬ 
notes the energy of the electron after the tunneling event. 
At zero temperature the Fermi distribution / becomes a 
step function and equations © and © have a simple 
interpretation: the electrons in the lead have no way to 
gain energy without absorbing a photon from the bosonic 
bath if fi = 0. This is why the convolution with the hot 
impedance is important - it permits the P(E) function 
to have some probability in the negative energy range, so 
a photon can be absorbed from that bath. It may then 
be given back to a cold bath as it tunnels left or right, 
as can be seen in Fig. 1(c). The asymmetry between 
the capacitances breaks the left-right symmetry. Equiva¬ 
lently we repeat that one may have taken the same value 
for Cl and Cr but with different impedances Zr and 
Zr. In this case of Fig. 1(b) the level is dephased by 
a fluctuating phase 0^, so that we assume a sequential 
non-resonant tunneling case and the rectified current is 
given bj®, 

/= T +>L (E d )T +>R (E d ) - r_, L (£ rf )r-, fi (£ d ) 
T+,l{-E d) + T +)R {E d ) + T_ )L (E d ) + T_ )R (E d ) 

( 7 ) 


B. Hot electrons and Zd — 0 (Fig. 1(b)) 

In this case no photons can be absorbed during the 
tunneling processes. However, we now consider the elec¬ 
trons to be thermally distributed in energy with a Fermi 
function so they can tunnel from lead to dot and vice 
versa. This symmetry can be broken by permitting the 
tunneling electrons to lose energy by emitting photons 
in either the left or right electrical circuits both tunnel¬ 
ing on and off the dot, emitting two photons, as can 
be seen in Fig. 1(d). Here, the asymmetry between the 
cold impedances breaks the left-right symmetry leading 
to rectification of current. 

As before we consider the case where the zero fre¬ 
quency external impedances of the transmission lines 
Z[ = Rj is small compared to the resistance quantum, 
Ri <C R q = h/e 2 , and define the large on = R q /Ri 1. 
In this case, we have again for long time, Ki(t ) = 
— (2/ai)ln(aiE c t/nh) + in/2 + q e . The rectified cur¬ 
rent assuming a non-resonant sequential tunneling is then 
given by Eq. 0, with the following expression for the 
tunneling rates, 


( nk B T V 

\qllE c J 


Li 


(- 


Tn±v/2 \ 

> k B T ) 


±n±n / 2 


T±, L (fi) = — To,z,e 

m ^ _2re (irksT\ a R T ( 

T ± , R m = -T 0 , R e (^r) ( 

where Li x (z) is the Polylogarithm function. Q denotes 
again the energy of the electron after the tunneling event. 


—e 




)■ 



FIG. 2. (color online). Rectified current with respect to 
Ea/ksT at zero bias. The red dotted curve corresponds to the 
setup displayed in Fig. 1(a), the dashed blue and the full blue 
curves correspond respectively to the non-resonant and the 
resonant case of the setup displayed in Fig. 1(b). We have 
ksT = 0.125 E c , oll(oll) — 5, olr{olr) — 30, To = 0.125E C 
and A;/r 0 = 1. Inset: Evolution of the maximum value of 
the current Im at the resonance with respect to Ec/ksT. 

However, in the case of Fig. 1(b), it is probably more 
relevant to address resonant tunneling processes. The 
transmitted electron can coherently bounce back and 
forth between the two barriers and interfere construc¬ 
tively before it exits. The rectified current is now / = 
(e/h) (r+ — T - ), obtained by summing over all the pos¬ 
sible paths of the electron, so that we have 

(r^+rM) 2 + (fi -^ )2 ’ 

r_, L (ft)r_, fl (n) 

( r°-qr °.« ) 2 + ( o —E d )2 

We have introduced r±j(Q) = HT±j(Q), and we will 
focus on the symmetric case where Tq,/ = = Tq 

in the following. This case permits Fabry-Perot type 
resonances between the two junctions forming the 
quantum dot 47 . It corresponds to a resonant sequential 
tunneling associated to the emission of two photons. 

Putting these results together, we present an analysis 
of the evolution of the rectified current with respect 
to the dot energy level Ed and to the charging energy 
E c . Here we fix the temperature T and To such that 
k B T/T 0 = 1. Note that T is the temperature of the hot 
impedance for the setup displayed in Fig. 1(a), while it 
is the temperature of the leads for the setup displayed 
in Fig. 1(b). The charging energy is the highest energy 
scale of the problem, and we do not calculate the current 
at temperatures of the order or greater than the charging 
energy. Above this threshold, the description we made 
is no longer accurate as the other levels of the quantum 
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FIG. 3. (color online) Left panel: Rectified current with re¬ 
spect to the left and right cavity frequencies ujl and ujr. We 
have k B T/r 0 = 1.2, Ed/To = 2, Z B lR q — ZR/R q — 0.5 and 
li — 0. Right panel: Rectified current with respect to ujr , for 
a setup with only one cavity coupled to one of the leads. 

dot must be taken into account. 


In Fig. 2, we plot the rectified current with respect to 
the dot energy for 0 ^( 0 ^) = 5 and aR(otR) = 30, at a 
fixed value of the charging energy E c /Tq = 8. We notice 
that this current is maximal at a dot energy of the order 
of the temperature. The value of this maximal current 
also depends on the charging energy of the dot. The in¬ 
set shows then the evolution of the maximum value of the 
current with respect to Ec/ksT (for E c /kRT > 2). We 
notice that the current decreases with the charging en¬ 
ergy, as expected from the power law evolution of P(E). 
The evolution of the current with respect to the dot en¬ 
ergy looks different for the case of Fig. la and the case 
of Fig. lb. These differences come from the effect of 
temperature, which contributes through a Bose-Einstein 
distribution for the case of Fig. la and through a Fermi 
distribution for the case of Fig. lb. The latter is respon¬ 
sible for the sharpness of the resonance (dashed and full 
blue lines) while the first one leads to a rather long tail 
(dotted red line). 


C. LC circuit and cavities 

Finally we focus on the setup displayed in Fig. 1(b) 
where the lead l is now coupled to a zero-temperature 
resonator of frequency uji = 1 /^LjCf (composed of an 
inductance Li and another capacitance Cf). We then 
have 


W) = Y q ( e ~ iMlt ~4 • 


( 10 ) 
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huj R /r 0 
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FIG. 4. (color online) Rectified current with respect to the 
temperature, and the cavity frequency. We have Fo,l = 
To ,r = To, Ed/To — 6, ZR/R q — 0.5. 


e K ‘ m =E~ 0 ELolr(;)(-l)"- 1 (f )V~"(11) 

This allows us to reach the following expression for Pi ( E ), 


Pi(E) = E^o h (f)* e“£ S(E - khun). (12) 

The left and right tunneling rates given by the Eqs. 
(5) and (6) read 


T ± , L (tt)=T 0 , L J2 

k 

T±, R m = T 0 ,RE 


(so 1 


kl 


e R i f(khu>L ± e) (13) 


(ft)' 


k\ 


e R i f(khcj R Tty- ( 14) 


These rates take into account tunneling events with the 
emission of several photons, and the ratio Zi/R q deter¬ 
mines the dominant processes occuring in the device. We 
first consider small values of Zi/R q where one-photon 
processes are the most relevant. 

Here we first fix the temperature T and To such 
that ksT/ Tq = 1.2. In Fig. 3 we plot the rectified 
current for different values of c jl and ujr , but with 
constant Z B = Zr. We have then photon-assisted 
tunneling and current rectification without any drive 
on the cavity, in contrast to Refs. I48H501 The best 
situation corresponds to very asymmetric configurations 
(blue and red regions), where one of the two cavities is 
suppressed. In this case, the coupling to a frequency- 
selective cavity leads to a rectified current substantially 
greater than in the case of two standard resistances in 
a highly asymmetric configuration (see Fig. 2 and Fig. 3). 


where we have again Z\ — y Li/Cf . We find 
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FIG. 5. (color online). Rectified current with respect to 
an external bias at different lead temperatures, ranging lin¬ 
early from the limit of zero temperature (full cyan curve) to 
ksT/To = 3 (full red curve). We have taken the resonant 
case of the setup displayed in the Fig. 1(b) (and adjusted the 
value of Ed/ksT at small bias). Inset: generated power with 
respect to the external bias. We have oll — 5, olr — 30 and 
r 0 = 0.125F C . 

We now study the dependence of the rectified current 
with respect to the temperature and the cavity frequency 
in this highly asymmetric case, where only the right lead 
is coupled to a cavity. In this case, the current is simply 
given by 


I = 


hR a 


J 


dPt 


F 2 
1 o 


n + (n - E d y 


F(Sl), (15) 


where 

F(ii) = [f (fi) f (hu R - n) - / (-0) / + fi)]. 

(16) 

The rectified current is written in Eq. ( p~5[ ) as a 
difference of two terms, which are products of two Fermi 
functions. At zero temperature both terms are zero 
(because one of the Fermi function is zero in both terms), 
and no tunneling events can occur. In the limit of high 
temperature, the flattening of the Fermi functions erase 
the asymmetry of the setup introduced by cjr, which 
also leads to a zero rectified current. We deduce the 
existence of one resonance: there exists an optimal lead 
temperature for which current is maximal. 

In Fig. 4 we plot the rectified current with respect to 
the temperature and to the cavity frequency. In this case 
of one-photon processes, we naturally find that the reso¬ 
nance occurs at a temperature T such that ksT ~ Jtljr. 
For greater values of Zr/R q , this resonance is shifted to 
a greater value of the temperature T, which can be esti¬ 
mated by ksT ~ khujR, where k is the number of photons 
of the dominant process at this value of Zr/R q . 



FIG. 6. (color online) Maximal power Pm with respect to 
ksT/To. The red curve with down-pointing triangles and the 
blue curve with circles correspond respectively to the resonant 
and the non-resonant case of the setup displayed in Fig. 1(b). 
The black curve with up-pointing triangles correspond to the 
setup displayed in Fig. 1(a) with Ai/To = 1. For these setups, 
we have F c /Fo = 20, oll [oll) — 5 and Oi?(d#) = 30, and Ed 
is adjusted to maximize the power. The green curve with 
squares corresponds to the case of one cavity, and ujr and Ed 
are adjusted to maximize the power. 

IV. CHARACTERISTICS OF HEAT ENGINES 



k B T/ r 0 


FIG. 7. (color online) Efficiency of the engines with respect 
to ksT/ ro- We use the same color and symbol legend as in 
Fig. 6. Parameters are unchanged. 


The generated power is given by the product I x V. 
When I > 0 and fi < 0, the device produces power: the 
generated current flows against a load potential. In Fig. 5 
we plot the evolution of the resonant current with respect 
to an external bias field at different lead temperatures, 
for the setup displayed in Fig. 1(b). We see that the 
increase of the temperature has two effects: it leads to a 
drop in the conductance, and to the appearance of the 
zero-voltage rectification current. 
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The heat engines can be characterized thanks to two 
physical quantities which are the maximum power that 
can be generated, whose evolution is shown in the inset of 
Fig. 5, and the efficiency of the engine. We plot in Fig. 6 
the maximum power as a function of temperature. In the 
single cavity case it exhibits a peak at a given tempera¬ 
ture, while on the other hand the power generated by the 
circuit coupled to the impedances is increasing with the 
temperature. In all cases we must consider thermal en¬ 
ergies which are smaller than the level spacing of the dot. 
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In the case of hot leads and cold environment, the effi¬ 
ciency is given by 77 = IV/(IV — J), where — J is the heat 
expelled to the cold environment, which is then given by 
the sum of all the energy carried by the emitted photons. 
The heat currents associated to the tunneling events are 


J±, L (n)=r 0 ,L J def(e±ttT fJ,/2)P L (e)e, (17) 

J±M n ) = r o ,R J def(e T ^ T n/2)P R (e)e. (18) 

Defining J(f2) = J +jL (f2) + J_ iL (0) + J +jL (fi) + 
J+, #(U), the emitted heat current is simply given by 
J(Ed) in the case of non-resonant tunneling, while it is 


J = [ dn J(S2)- 

J fr o 


^ IMl+To^^ 


( r °-q r °- fl ) +(n-E d )i 


(19) 


in the case of resonant tunneling. On the other hand, for 
cold leads and a hot environment, the efficiency reads r] = 
IV/ Jh, where Jh is the amount of heat received from 
the environment. We plot in Fig. 7 the efficiency of the 
Nano-engines with respect to ksT/To. We can remark 
that the efficiency of the setup displayed in Fig. 1(b) is 
similar for the non-resonant and the resonant case. 


Appendix A: Voltage fluctuations and the 
environment 


v L 


C c L : 


Cl 
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FIG. 8. (color online) Left impedance in the case of Fig. lb. 
The Left impedance is here a collection of capacitors and in¬ 
ductors. The correspondence between impedance and trans¬ 
mission line just fixes the ratio between inductances and ca¬ 
pacitances. We choose the capacitances to be equal to (7£ and 
adjust the inductances such that the ratio \JLl /C£ gives the 
resistance at uj = 0. 


V. CONCLUSION 

We have shown that the quantum vacuum fluctua¬ 
tions of a T = 0 electrical circuit permit the rectification 
of electrical current in a mesoscopic quantum system. 
While external energy is required for this rectification to 
happen, the fact that the electrical environment absorb 
photons and electrons differently for different external 
impedances or capacitances in the system is sufficient 
to break the left/right spatial symmetry and rectify the 
current. Of course, there are other ways to break this 
symmetry, such as having a nonlinear electron transmis¬ 
sion that is asymmetric (different tunneling rates is not 
enough). However, we stress that for our system in terms 
of the electron degrees of freedom, the system is com¬ 
pletely left/right symmetric - it is the vacuum fluctua¬ 
tions from the bosonic degrees of freedom that break the 
symmetry. 


We consider the case pictured in Fig. 8 (coupling be¬ 
tween the Left impedance and the voltage fluctuations 
on the left capacitor), where a resistance is modeled by 
an infinitely long transmission line. The Hamiltonian of 
the environment readJ^ 



where n denotes the site index of the transmission line. 
Following Ref. [51] we can diagonalize the Hamiltonian by 
putting 


Qn 



dx cos 

dx cos 


^ 7 TUX + 
( 7TTIX + 


f )om 
”)*>• 


(A 2 ) 



























We can verify that the couple of operators Q{pc) and c/)(x) 
satisfy the canonical commutation relations. Introducing 
then creation and annihilation operators aj, and a x 


4>fa) =7 


i^/2Cffaw a 

A( \ I C l(* 
Q 0 ) o — fa 


fax tile) 5 

+ <*x) > 


(A3) 


and the tunneling element from the left lead to the dot is 
dressed with the phase factor </>£ (in the main text) which 
corresponds to — 0o using the notations of the Appendix. 
We have for the time dependent correlation function of 
the variable </>o 


K L (t) = 


f 


e 2 cos 2 (7tx/2) ( _• t 

dx - fac fa 

huj x C c L 


- 1) • (AS) 


with uj x = 2 sin(7rx/2)u; c and uj c = I/^/LlC^ we get 



The voltage fluctuations on the junction capacitor are 
coupled to the environmental degrees of freedom and we 
have eSVL = eQo/C £ = f dx\ x (aj. + <Uc), where = 
cos(7rx/2). This coupling term is absorbed by 
the unitary transformation introduced in the main text, 


After a change of integration variable from x to w x , we 
finally reach 



(A6) 


The behaviour of this function is dominated by the low- 
frequency modes uj x /uj c <C 1, so that we can consider 
equations (3) to be valid with Zl = y/LiJ • 
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